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0.  Preliminaries:  Suppose  a  lot  of  transistors  is  to  be  divided  among 
two  purchasers.  The  lot  is  dichotomized  by  some  quality  criterion  and 
let  us  assume  there  are  m  high  quality  items  and  n  other  items  of 
tolerable  quality.  Purchaser  G,  who  has  contracted  for  r  items, 
suspects  the  supplier  of  favoring  purchaser  H  with  respect  to  the  m 
high  quality  items.  He  therefore  decides  to  test  the  hypothesis  that 
the  lot  has  been  divided  in  a  random  manner  against  the  alternative  that 
the  supplier  is  biased  in  favor  of  purchaser  H.  Fisher's  exact  test 
is  appropriate  to  test  the  null  hypothesis  but  in  order  to  examine  the 
sensitivity  of  the  test  one  must  define  and  parametrize  "degree  of 
bias"  and  obtain  the  non  null  distribution. 

This  report  gives  an  intuitively  appealing  meaning  to  "degree  of 
bias"  when  sampling  from  a  finite  population  and  obtains  a  distribution 
herein  called  the  noncentral  hypergeometric  probability  distribution. 
Various  properties  of  the  distribution  are  studied  and  asymptotic  results 


given. 


SECTION  1 


INTRODUCTION 

1.1.  The  problem:  Consider  a  set  S  of  i  elements  which  is  dichotomised 
in  some  manner  (say,  by  an  observable  characteristic)  into  subsets  M 
and  N  containing  m  and  n  ■  A  -  m  elements,  respectively,  and  a 
sampling  mechanism  vhich  in  some  way  selects  a  subset  R  of  r  elements 
from  S.  Let  a  be  a  realization  of  the  random  variable  A  denoting 
the  number  of  elements  in  ROM.  The  results  of  such  an  experiment  are 
usually  displayed  in  the  form  of  a  2x2  contingency  table  such  as 


Table  1.1 


a 

c 

m 

b 

d 

n 

r 

s 

l 

where  b  =  r  -  a 

c  =  m  -  a 
d  =  n  -  r  +  a 
s  =  A  -  r 

We  wish  to  test  the  hypothesis  that  the  sampling  mechanism  is  random 
against  the  alternative  of  "non-randomness".  The  notion  of  random 
sampling  from  a  finite  population  is  well  defined  by 

Definition  1.1:  A  sample,  each  of  whose  elements  is  drawn  from  a  finite 
population,  i6  said  to  be  a  random  sample  if,  at  each  draw,  all  elements 
available  for  selection  have  an  equal  probability  of  being  selected  [8]. 


Under  the  hypothesis  of  randomness,  the  random  variable  A  has  a 


hypergeometric  distribution  with  parameter*  a, a  and  r  and  bene*  Fisher's 
“exact"  test  [6]  is  applicable,  i.e.  the  hypothesis  of  randomness  is 
accepted  if  and  only  if  L  <  a  <  U  where  L  and  U  are  percentage 
points  of  the  appropriate  hypergeooetric  distribution  and  depend  on  the 
specified  significance  level  a.  The  alternative  of  “non- random 
selection",  however,  is  not  well  defined  and  hence  discussions  of  the 
sensitivity  of  the  test  are  impossible. 

The  purpose  of  this  paper  is  to  define  non- random  or  biased  selec¬ 
tion  in  an  intuitively  appealing  way  and  to  parametrize  the  degree  of 
bias  in  order  to  obtain  the  non-null  distribution  of  the  test  statistic 
A  and  thereby  determine  operating  characteristics  of  Fisher's  “exact" 
test. 

In  order  to  relate  the  model  under  consideration  to  other  models 

which  yield  data  in  the  form  of  table  1.1,  a  brief  survey  of  statistical 
inference  based  on  2X2  tables  is  given  below. 

1.2.  2X2  tables:  Several  types  of  statistical  investigations  lead  to 
data  dlsplayable  in  the  form  of  table  1.1.  The  particular  "table" 
obtained  as  the  result  of  a  random  experiment  can  be  compared  to  the  set 
of  all  realizable  tables  and  thus  provide  a  basis  for  inference.  G.  A. 
Barnard  [1]  and  E.  S.  Pearson  [9]  have  suggested  classifying  models  by 
the  number  of  fixed  marginal  totals.  For  example,  if  m,  n,  and  r  are 
fixed  then  specifying  a  value  A  *  a  completely  determines  the  other 
entries  so  that,  in  this  case,  the  sample  space  is  one  dimensional.  If 
r  is  not  fixed,  then  the  sample  space  is  two  dimensional  since  one  needs 
to  specify  a  value  (A,B)  ■  (a,b)  to  completely  determine  the  table. 

These  classifications,  together  with  the  names  used  by  Barnard  and 
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Pearson  are  given  belov  in 


Table  1.2 


Fixed 
Parameters 


Classification 


Barnard 


Pearson 


Sample  Space 


m,n,r 


Independence 

trial 


n 

■ 

ca 

■ 

■ 

D 

o 

■ 

■i 

->  A 


max{0,  r-n)  <  A  <  min  (r,m) 


m.n 


Comparative 

trial 


II 


A 

m 

B 

n 

->  (A,B) 


0<A+B<m+  n>A  <  m 


Double 

Dichotomy 


III 


D 

■a 

■ 

m 

H 

■ 

■ 

D 

->  (A,B,C)| 


0<A  +  B  +  C<i 


Classically,  the  type  II  table  is  pertinent  to  the  problem  of  com¬ 
paring  binomial  parameters  of  2  large  populations  from  which  random 
samples  of  sizes  m  and  n  have  been  drawn  while  the  type  III  table  is 
used  in  connection  with  testing  a  single  large  population  for  indepen¬ 
dence  of  2  characteristics  based  on  a  sample  of  size  i.  In  these  two 
cases,  uniformly  most  powerful  unbiased  tests  exist  [7]  and  the 
sensitivity  of  such  tests  are  discussed,  for  example,  in  [2],  [4],  [11], 
[12]  and  [14] .  Type  I  tables  are  frequently  employed  in  connection  with 
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tasting  ths  hypothesis  that  two  treat— nts,  say  sad  T^,  in  iden¬ 
tical  in  so  far  as  producing  *  "reaction"  In  a  croup  of  I  individuals 
which  has  been  randosly  divided  into  two  groups  so  that  ■  raoeive 
and  n  «  i  -  a  receive  •  Without  further  assusqptions,  the  inferenoee 
pertain  only  to  the  reactions  of  the  group  being  tested.  Indeed,  it  way 
be  impossible  to  repeat  the  experiment  if,  for  example,  reaction  Inplles 
death  or  some  other  less  abrupt  change  in  susceptibility.  Pearson  [10] 
considers  a  particular  one-sided  sensitivity  analysis  which,  however. 

Is  not  applicable  to  the  problem  under  consideration. 


SECTION  2 


BIASED  SAMPLING 


2.1.  Definition  and  Parametrizatlon:  Imagine  a  sampling  mechanism 
sequentially  drawing  a  sample  of  size  r  from  the  set  described  in 
section  1.1.  Prior  to  the  k  draw,  k=l,...,r,  let  m(k)  and  n(k) 
denote  the  number  of  elements  remaining  in  sets  M  and  N,  respectively, 
and  let 


X(k) 


i  if  tth 
0  otherwise 


draw  is  an  element  of 


M 


A  =  I  X(k)  . 
k=l 

Definition  2.1;  Hie  sampling  mechanism  is  biased  by  9,  0  <  6  <  «, 

if  and  only  if 


(2.1.1) 


Remarks: 

(1)  Replacing  P(X(k)  =0}  by  l-P{X(k)  =1}  we  obtain 

?«<*>  ■  u  -  .w  .  (2-1-2) 

(2)  If  m(k)  =  n(k)  then  P(X(k)  =  0)  =  9  P{X(k)  =  1)  so  that 
definition  2.1  has  the  intuitively  desirable  property  that  when  the 
remaining  sets  are  equal  in  size,  a  bias  of  6  implies  the  sampling 
mechanism  is  "9  times  more  likely"  to  select  an  element  from  N  than 
from  M. 
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(3)  If,  «t  the  kth  draw,  m(k)  ■  0  n(k)  then 
P(X(k)  ■  1}  -  P(X(k)  >  0}  -  1/2.  Bence,  even  though  the  sampling 
mechanism  is  biased  by  0,  it  is  equally  likely  that  the  draw 

will  yield  an  element  from  M  as  an  element  from  R  when  the  elements 
remaining  in  M  are  0  times  more  numerous  than  those  remaining  in  R. 


(4)  lim  P(X(k)  -  1} 
9-*  “ 


)_{° 

th  , 


if  m(k)  *  0 
if  m(k)  /  0 


Hence,  for  a  large 


degree  of  bias,  the  k  “  draw  will  be  an  element  from  set  M  provided 
M  has  not  been  totally  depleted. 

(5)  For  0  =  1,  definitions  2.1  and  1.1  coincide  and  clearly  A  has 
a  hypergecmetric  probability  distribution  with  parameters  m,  n,  and  r. 
For  this  reason,  the  distribution  of  A  for  0  <  6  <  ®  will  be  called 
noncentral  hypergeometric  with  noncentrality  parameter  9  and  will  be 


denoted  P0{A  =  a). 


2.2  Example:  Hie  main  result  of  this  paper  was  motivated  by  the 
following  genetics  study  [ 3] .  A  small  rocky  island  located  about  30 
miles  seaward  of  San  Francisco  Bay  has  supported  a  population  of  rabbits 
for  at  least  100  years.  This  population  affords  the  geneticist  an  ideal 
nonlaboratory  environment  in  which  to  study  genetic  fixation  ^  because 


1/  Genetic  fixation  in  a  specified  population  is  the  tendency 
for  all  but  one  allele  to  disappear  from  a  particular  gene 
locus.  [13] 
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(a)  Migration  of  animals  in  nonexistent 

(b)  Die  population  is  greatly  reduced  each 
Sumner  by  the  sparce  supply  of  food  ^ 

(c)  Absence  of  predators 

(d)  Hie  terrian  facilitates  trapping  to  the 
extent  that  it  is  safe  to  assume  that 
sampling  error  is  negligible. 


One  of  the  loci  under  consideration  is  blood  type.  If  no  genetic 
selection^has  taken  place  in  the  population,  the  number  of  generations 
and  the  effective  size  of  the  breeding  population  can  be  used  to  deter¬ 
mine  the  probability  of  fixation  [5].  For  the  case  in  question,  this 

probability  is  approximately  0.04.  Since  the  blood  type  locus  was 

4/ 

observed  to  be  polymorphic,-'  a  study  was  made  tc  determine  if  selection 
was  occuring  in  neonate  survived.  Neonates  were  captured,  blood  typed, 
released  and,  after  three  months,  the  surviving  animals  were  recaptured. 
The  data  are  given  below  in 


Table  2.2 


Survived 

Died 

Totals 

Homozygotes 

41 

5^ 

95 

Heterozygotes 

54 

86 

120 

Totals 

75 

140 

215 

2/  In  1961,  for  example,  only  69  of  1200  neonate  rabbits  sur¬ 
vived  to  an  age  where  reproduction  could  occur. 

2/  Changes  in  allele  frequence  between  generations,  other  than 
those  attributable  to  random  gamete  and  survival  sampling, 
are  called  selective  changes.  [15] 

kj  More  than  one  allele  present  in  the  population  at  the  locus 
in  question;  the  opposite  of  genetic  fixation.  [13] 
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The  hypothesis  to  be  tested  is  that  the  sample  of  size  75  was  drawn 
randomly  by  "nature"  from  the  population  of  size  215  against  the 
alternative  that  nature  is  biased  (i.e.,  selection  against  one  phenotype 
or  the  other  by  the  mechanism  of  death  had  oc cured) .  To  test  this 
hypothesis,  one  simply  applies  Fisher's  exact  test  at  an  appropriate 
level  using  the  values 

l  m  =  95 
n  =  120 
r  =  75 
A  *  41 

« 

Attempting  to  parametrize  the  alternative  and  obtain  the  power  function 
of  the  test  motivated  this  paper.  After  the  noncentral  hypergecanetric 
distribution  is  obtained  in  section  5,  it  will  be  applied  in  section  5 
to  a  discussion  of  the  power  of  the  above  test. 
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SECTION  3 


THE  NONCENTRAL  HYPERGEOMETRIC  DISTRIBUTION 


3.1.  A  First  Expression;  In  this  subsection  we  shall  obtain  an  expres¬ 
sion  for  Pfl{A  ■  a}  by  a  straight  forward  enumeration  method.  This 
form,  although  not  mathematically  appealing,  will  be  used  in  a  discussion 
of  bounds  on  P.{A  =  a)  and  in  certain  asymptotic  results  in  section  4. 

9 

A  simpler  and  more  desirable  expression  will  be  obtained  in  subsection  3-2. 

Using  the  notation  introduced  in  Section  2,  let  X  =  (X(l), 

X(2),  ...,  X(r))  be  the  random  variable  denoting  the  outcome  of  drawing 
r  elements  from  MON  and  x  =  (x(l),  x(2),  ...,  x(r)}  be  a  realization 
on  X.  Then,  using  (2.1.2)  and  P{X  =  x)  *  P{X(l)  *  x(l)}  x 
P(X(2)  =  x(2)  |  X(l)]  P(X(3)  =  X(3) j  X(l) ,  X(2) )  ....... 

P(X(r)  =  x(r) |  X(l),  X(2),  -  X(r-l))  we  obtain 


/_£ 


]  |  x(k)  m(k,x)  +  [l-x(k)]  8  n(k,x) 
m(k,x)  +  8  n(k,x) 


k=l 


P{X  =  x]  =< 


max  (o,r-n)  <  £  x(k)  <  min  (r>m) 
“  k=l 


otherwise 


(3-1.1) 


where  m(k,x)  and  n(k,x)  are  the  number  of  elements  remaining  in  M 
and  N  respectively  prior  to  the  ktlx  draw.  Actually,  m(k,x)  and 
n(k,x)  depend  only  on  x(l),  x(2),  ...,  x(k-l)}  more  precisely, 

k-i 

m(k,x)  =  m  -  Y  X(J) 

J-l  (3.1.2) 

k-i 

n(k,x)  =  n-  7  [l-x(j)]  , 

J-l 
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Then,  P.(A  ■  a)  ■  £  p(X  ■  x}  where  X(a)  =  {x 

xc  X(a) 


l  x(k)  -  a)  (3.1.3) 
k«l 

Viewing  the  (*)  elements  in  X(a)  as  binary  numbers  instead  of 
vectors  for  a  moment,  let  us  order  X(a)  by  decreasing  magnitude  and 
index  the  ordered  set  with  J . 

i.e.  J  *  1  (l#  1>  •••  1,  1,  0,  0,  ...  ,0) 

J  *  2  (l,  1,  ...  1,  01,  0,  ,  0) 


J  -0 


(o,  0,  ...  0,  1,  1,  •••  ,1) 


where  there  are  exactly  a  l's  and  r-a  0's.  Denote  the  j  x  in 
the  ordering  of  X(a)  by  Xj .  By  using  (3.1.1)  and  (3.1.3)  and  obser¬ 
ving  that  for  each  x  c  X(a), 

J 


T T  [xj(k)  m((k>xj)  +  [1  -  Xj]  6  n  (k,  Xj)J  =  (m)a(n)r  a9r‘a  where 


(m)  =  m(m-l)  ....  (m-a+1)  we  obtain 

& 


P0{A=a)  =  (m)a(n)r_aS 


r-a 


(?) 


a)  r  f 

I  TT  U(k,x  )  +  0n  (k,x  ) 
1=1  k=l  L  3  3 


-1 


(3*1.10 


As  a  verification,  note  that  for  0=1 


TT  £n(k,Xj)  +  0n  (k,x,)J  =  |(i)rJ 


(m)  (n) 

so  that  P{A=a)  =  &^) ' ~ 


which  is  the  central  hypergeometric 


case. 
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3-2  A  Preferred  Expression;  While  (3.1.4)  is  correct,  it  is  a 
complicated  notational  expression  and  leaves  much  to  be  desired.  The 
attack  used  in  this  subsection  is  that  of  partial  difference  equations. 
In  order  to  write  down  a  partial  difference  equation  we  shall  need  to 
modify  the  notation  for  Pg(A=a)  in  order  to  include  reference  to  the 
particular  values  of  the  parameters .  Therefore,  let  us  replace 

P0(A=*a}  by  P0{a;  ra,  n,  r). 

Conditioning  on  X(l)  we  obtain 


P0{a;  m,  n,  r)  =  E  P0(a;  m,  n,  r 


*  PQ{a;  m,  n,  r 


X(l))j 

X(l)  =  1}  P{X(1)  =  1} 

X(l)  =  0)  P{X(1)  =  0} 


+  P0(a;  m,  n,  r 


Pe(a-1;  m-1,  n,  r-1]^  +  P0(a;  m,  n-1,  r-lJjgL 


Multiplying  by  m  +  n0  we  obtain 
(m+nQ)  P0(a;  m,  n,  r}  =  m  PQ{a 

with  boundary  conditions 


m-1, 

n. 

r-l)  +  nS  P0(a;  m. 

n-1,  r-l} 

(3.2.1) 

rH 

'0 

.1 

a  >  0 
a  =  0 

(3.2.2) 

r]  = 

ro 

[1 

a  <  r 

a  =  r 

(3-2.3) 

0}  =  ■ 

& 

a  >  0 
a  =  0 

(3.2.4) 

r}  = 

( 

n) 

r 

(3.2.5) 

(2 

+  n)r  • 
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The  first  three  boundary  conditions  need  no  explanation*  the  fourth  is 
obtained  directly  from  (3-1.4). 

Now  (3.2.1)  is  a  partial  difference  equation  in  four  variables. 
Often,  by  assuming  the  discrete  variables  are  continuous  and  viewing  the 
partial  difference  equation  as  a  partial  differential  equation,  a  simpli¬ 
fication  of  the  former  will  be  suggested  by  the  form  of  the  solution  of 
the  latter.  Rewriting  (3.2.1)  as 

m(Pg{a;  m,  n,  r)  -  P0{a-lj  m-1,  n,  r-1))  +  n©  (P0(a*  m,  n,  r)  - 

P0iaj  m,  n-1  r-1))  =  0  and,  by  adding  and  subtracting  appropriate 
terms  we  obtain 


m(PQ{a;  m,  n,  r)  -  P^aj  m-1,  n,  r )) 

+  m(PQ{a;  m-1,  n,  r)  -  P0(a;  m-1,  n,  r-1)) 

+  m(P0{a;  m-1,  n,  r-1)  -  P0[a-lj  m-1,  n,  r-1)) 
+  n©  (P{aj  m,  n,  r)  -  P0{a>  m,  n-1,  r)) 

+  n©(p{aj  m,  n-1,  r)  -  P(a*  m,  n-1,  r-1))  =  0 
which  leads  to 


Applying  the  method  of  characteristics, 

dm  _  dr  _  da  _  _dn  _  dp 
m  -  m+  n©  m  "  n©  0 

The  first  and  third  members  of  (3*2.8)  yield 


(3.2.6) 


(5.2.7) 


(3-2.8) 


m 


a  +  c 


1  . 


(3.2.9) 
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The  first  and  fourth  members  of 

In  n  * 
or  n  = 

The  first  and  second  members  of 

dr 

dm 


r 

Hence,  the  solution  of  (3.2.7) 
or,  using  (3-2.9),  (3-2.10)  and  (3-2.11) 

P  {aj  m,  n,  r]  =  F(m-a,  nm  ,  r-m-n)  ,  (3-2.12) 

U 

This  result  suggests  the  following  change  of  variables.  Let 

P  la;  m,  n,  r)  =  Q(m+n-r,  m-a,  m,  n) .  (3-2.13) 

& 

Then  (3-2.1)  becomes 

(m+n0)  Q  (m+n-r,  m-a,m,n)  =  m  Q  (m+n-r,  m-a,  m-1,  n) 

+  n0Q  (m+n-r,  m-a,  m,  n-l)  ,  (3-2. 14) 

Notice  that  the  first  two  arguments,  (m+n-r)  and  (m-a)  are  Identical 
in  each  term  of  (3-2.14)  so  that  the  number  of  variables  has  been  re¬ 
duced  from  four  to  two.  In  terms  of  Q,  the  boundary  conditions 


(3-2.8)  yield 


0  In  m  +  In  C0 


(3-2.10) 


c2m 


0 


(3-2.8)  yield 


m+n0 

m 

1  + 

n0 

m 

1  + 

0C2m0_1 

m  + 

C2m0  +  C 

m  + 

n  +  0^  . 

(3-2-11) 

will  be  some  function  of  (c^,  c2,  c^) 
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(3.2.2)  through  (3*2.5)  become 

Q(n-r,  -  a,  0,  n)  -  [®  J  ^  ® 

Q(°-ri  m-a,  m,  0)  -  ' 

Q(m+n,  m-a,  m,  n)  =  jj  *  >  ° 

(n) 

Q(m+n-r,  m,  m,  n)  = 


(5  +  »>, 


(3*2.15) 

(3-2.16) 

(3.2.17) 

(3.2.18) 


Now,  define 


m+n-1  min^r.m) 
r=0  a=max(o,r-n) 


r  a 


®m  n(x»y)  =  I  1  Q(m+n-r,  m-a,  m,n)x*  y 

J  %• _ i 


(3.2.19) 


Boundary  conditions  on  4>  (x,y)  are  easily  obtained  from  (3-2.15)  and 

iQy  n 

(3.2.16). 


n  , 

%  _(x,y)  =  l  Q(n-r,  0,  0,  n)  xr  =  £-=±  (3.2.20) 

u>  r=0  x  " 

»  Q(x,y)  =  I  Q(m-r,  m-r,  m,0)  xr  yr  =  (3.2.21) 


Multiplying  (3-2.19)  by  m+ne  and  splitting  off  the  term  for  r=0  we  obtain 

min{r,m) 

I  < 

a=max(o,r-n) 


m+n-1  min{r,m) 

(m+ne)  <J>  (x,y)  =  V  T  (m+ne.)  Q  (m+n-r,m-a,m,n)x  y  +  (m+ne) 

m, n  ' 


and,  applying  (3-2.14)  to  the  right  hand  side,  we  obtain 

n-1  min(r,m) 

i,  tl  , 

=1  a=maxlo,r-nJ 


m+n-1  min(r,m) 
(m+ne)  n(x,y)=  m  £ 

;  r= 


X*  & 

2_  Q(m+n-r,  m-a,  m-1,  n)x  y 


m+n-1  min{r,m) 

+  n0  Y  Y.  Q(m+n-r,  m-a,  m,  n-l)xr  ya 

r=l  a=max{o,r-n} 


+  (m  +  n6) 


(3.2.22) 
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In  the  first  term  of  the  right  hand  side  of  (3.2.22),  max{o,  r-n)  can 
be  replaced  by  max{l,  r-n}  since  Q(m+n-r,  m,  m-l,n)  *  0.  Then,  re¬ 
placing  a  by  a-i-1,  we  obtain 
m+n-1  min{r-l,m-l} 

“X  Z  Q(m-l+n-[r-l],  m-l-a,  m-l,n)xr  ya+ 

r=l  a=maxlb,r-l-n} 

and,  replacing  r  by  r+1,  we  have 

m-l+n-1  min{r,m-lj 

m  X  X  Q(m-l+n-r, m-l-a, m-l,n)xr+  ya+  =  mxy<X>  1  (x,y) 

r=0  a=maxl0,r-n}  m-i,n 

(3.2.23) 

In  the  second  term  of  the  right  hand  side  of  (3.2.22),  min{r,m}  can 
be  replaced  by  min(r-l,  m}  since  Q(m+n-r,  m-r,  m,  n-1)  =  0.  Then, 
replacing  r  by  r+1  we  obtain 

m+n-l-l  min|'r,m] 

n0  X  I  Q(m+n-l-r,m-a,m,n-l)xr+1  ya  =  n 0x<t  ,  (x,y) 

r=0  a=maxlO,r-[n-l] }  m,n-l 

(3.2.24) 

Substituting  (3-2.23)  and  (3.2.24)  into  (3.2.22)  yields 


(m+ne)  (x,y)  =  mxy  \_1>n  (x,y)  +  n  9  x  *  ml  (x,y)  +  (m+n0)  , 

(3.2.25) 

Let  \|r(|,T])  be  the  generating  function  of  4>  (x,y): 

m,  n 


m  n 


MS*1!)  =  X  X  ®  (3.2.26) 

«  «  m»n  m!  nl  •  '  7 

m=0  m=0  7 


Boundary  conditions  on  t(l#Tj)  are  easily  obtained  from  (3.2.20)  and 
(3.2.21)  as 

♦(cm)  -  £  v  £ 
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t  -  /  n  ,  v  n 

1  r-  ix_ 

*-1  nto  ni 


*X1 .  »1 
x-1 


«•*")  ■  z  Vo  S 


m=0 


(3.2.27) 


5=1  l  ('">■  -  1 


e^j-  e1 

xy-1 


(3-2.28) 


From  (3*2.26)  one  obtains  by  differentiation: 


00  00 


n 

1_!L 


i  ^ .  ( t . ti)  =  y  y  m  4  * ,  ^ , 

s  ’V5'  17  L  e-  m.n  min! 

b  m=l  n=0 


GO  00 


t)©  ^  (1,11)  =  y  y  n©  $  v  *, 

1  w  £.  i.  m  n  mini 

1  m=0  n=l  ’ 


m  n 

-2—3-  and  hence 


„m  n 

Ll 


!  t.(!>n)  +  n  0  *.(6^)  =  I  l  (ra+n0)  ®  n  Jnrr 

&  1  m=l  n=l  ’ 


.m  00 


y  m  $  +  y  n©  $  (3*2.29) 

mtl  ”<°  °'n  n: 


Substituting  (3.2.25)  into  the  first  term  on  the  right  hand  side  of 
(3.2.29)  yields 


.m  n 


.ran  00  00 


m  n 


xy  y  y  m$  .  +  x©  y  y  n<J>  .  +  y  y  m  Jj-A- 

m=l  n=l  m_1'n  mln-  m*l  n=l  m'n"1  m’n-  m=l.  .ntl  m'n’ 

“  “  sV  (3.2.30) 

+  ©  y  y  n  -v-v 

min! 

m=l  n=l 
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Taking  the  terms  of  (3.2. 30)  one  at  a  time: 


00  00 


m  *m-l,n  mini 

m=i  n*l 


.m+1  n 


xy  y  i  •  »  ,  -? 

^  in. n  ml  ni 

m=0  n*l  7 


w  I  I  v.  £r2r  -  ■*««.<» 


m=0  n=0 


=  |xyi(S,n)  -  xy| 


(3.2.31) 


00  00 


xG  X  I  n<J>  . 

mtl  n=l  m,n-1  mlnl 


oo  oo  m  n+1 

■  x0  £  7  0  _  ±rrr~ 

*-v  m.n  mini 
m=l  n=0  ’ 


xOr]  E  I  •  „  -  xerrtr(0,tj) 

m.n  mini 

m=0  n=0 


=  X0T|\|r(|,T])  -  X0T] 


(4^) 


(3-2.32) 


w  w  m  II 

l  T  m  VV 

S  mln! 

m=l  n=l 


00  00  m+1  -  n 

V  7  * - 3- 

S  ml  n! 

m=0  n=l 


=  {(e^  - .') 


(3.2.33) 


IT 


(3.2.5^) 


00 

®  E 

00 

I 

n 

n  JL2L 
mini 

m»l 

n*l 

00 

“  ,m 

L  ml 

= 

9 s, 

m»l 

n=0 

The  second  term  on  the  right  hand  side  of  (3-2.29),  after  applying 
(3-2.21)  becomes 


Similarly,  the  third  term  on  the  right  hand  side  of  (3.2.29),  after 
applying  (3-2.20),  becomes 


e 

x-l 


(3.2.36) 


Substituting  (3-2.31)  through  (3.2.3 6)  into  (3.2.29)  and 
simplifying  the  resulting  expression  yields? 


5  Vg(fiji))  -  xjnKl,?}) 


+  en  k^d/n)  -xt(!,n) 


=  (5  +  ©n)e5+n 
(3-2.37) 
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We  now  seek  a  solution  of  (3. 2. 37)  which  satisfies  the  boundary  condi¬ 
tions  (3.2.27)  and  (3.2.28). 


Let  UU,H) 


Then  (!,t|)  =  e_Xy5  '  XT*  ^(5>n)  -  xjn|r(|,T))j 

^  (i>T))  =  e-xy?  "  XT)  ^(£,1)  -  xijf( £, tj)J 

and  (3.2.37)  becomes 

£Ug(£>n)  +  ©^(£,1)  *  (£  +  ©ri)e"6(xy_1)"T,(x"1)  #  (3.2.38) 


The  boundary  conditions  become 


U(o,n) 


i.e-n^-1) 

x-1 


U(£,0) 


i.e-t(»y-1) 

xy-1 


(3.2.39) 


(3-2.40) 


The  Mellin  transform  of  a  function  U(£  ,r\)  with  respect  to 
( 5  will  be  treated  as  a  parameter)  is  defined  to  be 

lj((£>s)  =  f  n8"1U(£/n)dTi  where 

Jo 


s  is  a  complex  variable  suitably  restricted  to  insure  convergence,  if 
the  transform  exists  at  all.  Taking  the  Mellin  transform  of  both  sides 
of  (3.2.38)  yields 
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-  00  -00  .  . 

i  J  ns’1U|(*'Ti)dTi  +  e  j  tj8-1  hur,(5>Ti)|  dTi 


,  e-5(xy-l) 


Simplifying,  we  obtain 


r"  s-i 

['Jo  1 


t*-1*'1  a.,  * 


IU5(S,s)  -  es'u(s,s)  =  e's(xy'i)  fs(x-i)'sr(s)  ♦  0(x-i)‘(s+1)r(6+i)j 

-0S-1 

and,  multiplying  by  the  integrating  factor  §  gives 

4  [resU(s,»)j  -  r98-1  «-s(xy-1)  fs(x-i)"8r(s)  ,  0(x-i)-'s+1)r(s*i| 


and  hence, 


U(l>s)  -  t0S  JJ  t'98'1  e"T(xy_1)  |r(x^l)"8T(s)+0(x-l)-^S+1^r(s+l)  dr 


r(s+1h 


where  C  is  a  constant  of  integration  which  may  be  a  function  of  tj. 
Since  the  mathematics  is  simpler  and  a  solution  satisfying  the  boundary 
conditions  is  obtained  by  taking  C=0,  let  us  do  so;  also,  let  us  make 

T 

the  change  of  variable  t  =  y  . 


K.(5>s)  =  t”08”1  e"1^"1^  jtt(x-l)'sr(s)  +  ©(x-l)‘ 


(s+l) 


r(s+i) 


Inverting,  we  obtain 


U(i,n)  =  J  e'D^7"1^  +  Ti^x“1^t0]  (i  +  ©nt0-^ 


(5.2.41) 
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-  f  e"tl(X'1)te0rite'1dt 


—  f1 

x-1  / 

■Jo 


e ^(x-yt^at 


B  JL  -n(x-i)te 


(l-e  1^)  which  checks  with  (3*2.39) 


U(§,o)  =  f 


e-{(«y-l)te  4t 


1  -S(xy-l)t 
xy-1 


^3j(l-e  which  checks  with  (3*2.40), 


Using  (3*2 .41)  we  can  now  recover  \|r(|,T]) 

*(S,l)*  f1  ^((l-xyjwxy]  *  t(l-x)te* 

Jo 


=  I  Z  If  [(i-vJt+xyAtl-x)  t0+x]n[|+0T)t0" 
m=0  n=0  Jo 


i  ")  „m  n 

1]dttV17  (3*2.1 
J  I  min!  v 


From  this  expression,  the  coefficient  of 


m  n 

is  seen  to  be 

mini 


21 


Vn(x'y) "  ml  [(1-xy)t  +  xyf_1  [(i-*)te  +  x]ndt 


+  nfl 


J  [(l-xy)t  +xyjm  ^(l-x)t9  +  xj  te_1dt 

(5-2.43) 


From  the  definition  (3.2.19)  of  $  (x,y)  we  observe  that 

m,n 


min{r,m)  o  , 

n)ya  =  i 


a=max 


V  Q(m  +  n-r,  m-a,  m, 

:TO,r-n) 


•  Jx,y) 


m^n 


rl  Sxr 


x=0 


Considering  the  first  term  on  the  right  hand  side  of  (3*2. 43),  and 
applying  Leibniz's  rule  for  differentiation  and  interchanging  the  order 
of  differentiation  and  intergration, 


Jr ■  Yy  J  [(i-*y>* +  xy]  [(i-x)te  +  x]  dt 

ox  Jo 


*»0 


m  r1  I  Y  /r\  (l-xy)t  +  xy  m  1  dr  k  (l-x)t0  +  x 

r‘  1  k=0  '  dx  dxr'k 


x=q| 


dt 


■  *  r  1  (> 

'o  k=0  ' 


£}  (m-l)k(y[l-t])ktn-1-k(n)r.k(tl-te])r-kt8<n-r+k) 


dt 


m 


r 

V 

L 

k=0 


Mk  <n)r-k  r1 

L  ki  r-ki  J„ 


(i-t)t(i.te)r-ktm-k-1  *  9(,i-r+k)dt 
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■  l[w  ft)  f 


te}r-ktm-k-l  +  6(n-r+k)dt  yk 


(3.2  .44) 


Similarly,  the  Becond  term  on  the  right  hand  side  of  (3*2,43),  when 
differentiated,  becomes 

£  I7m\  /n-1^  f1  ,,  ^.k,,  ^r-k  m-k-1  +  6(n-r+k) 1  k 


£  [(-)  (£)£  (i-t)k(i-t9)r'k  t""*-1  +  j  A 


(3-2.45) 


Here,  it  is  understood  that  =  0  if  b>a.  By  combining  (3.2.44) 
and  (3.2.45)  and  by  uniqueness  of  power  series,  we  find  that 


(ajm,  n,r }  .  [  *  08  fa)  fr^]  X 


f1  +  8<”'r  *  S>4t 

Jo 


Some  manipulation  then  results  in: 


P0(ajm,n,r)  =^^|m-a+e(n-r+a)J  J  (l-t)a(l-t6) r  V1  a  1  ^  d^ 

(3.2.46) 

3.3  Compliments  and  verifications :  In  this  subsection  we  shall  verify 
that  (3-2.46)  has  required  properties  of  a  bonifide  probability  dis¬ 
tribution  and  satisfies  the  defining  conditions  (3.2.1)  through  (3-2.5) 

Theorem  3-3.1:  For  r  <  m+n  and  max(o,r-n)  <  a  <  min{r,m),  P0{ajm,n,r) 
as  given  in  (3.2.46)  is  discrete  probability  mass  function. 


Proof: 


Clearly  P  (aj  m,  n,  r)  >  0.  The  generating  function  for 


mlnCr, m} 

asmaJro,r-n) 


Pgfa;  m,  n,  r) 
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can  be  obtained  by  setting  y  ■  1  in  (3.2,19) 


m+n-1  I  min{r,m) 

Vn'*'1’  *  I  , 

r=0  |a=max 


I  ,  Va>  m,  n,  r))x* 

To/r-n)  |  (3*  3*1) 


=  m 


% 

f  1  .  ,m-l  r  -in 

J  [(l-x)t  +  xj  [(l-x)t0  +  xj  dt 


/■I,  lm  r  8  -in-i 

+  n  9  1  [(l-x)t  +  xj  [(l-x)t  +  xj  t  dt 

Jo 


(3-3.2) 

by  (3-2.43). 

Integrating  the  first  term  on  the  right  hand  side  of  (3-3.2)  by  parts: 
Let  U  =  [(l-x)t0  +  x]n,  dV  =  m[(l-x)t  +  x]m  ‘'’dt 

then  dU=n[(l-x)tS+x]n"10(l-x)t0"1dt,  V  =  [(l-x)t  +  x] 


we  obtain  the  expression 


Tji^y[(l-x)t0+x]  [( 1-x) t+x]  -  nS  f  [(l-x)t+x]  [(l-x)t0+x]  t0'1dt 

'  '  O  -Jo 


which,  when  substituted  in  (3 -3.2)  gives 


«  (x,l)  .  fi  =7 

m,nx  1-x 


,  m+n  m+n-1 
1-x  r  r 

x 


(3.3.3) 


Since  (3.3-1)  and  (3.3.3)  are  identities  in  x,  we  have 


min{r,m)  m  rl  - 

“  (“)(  ")[m-a+0(n-r+a)] 

a=maxTo,r-n]  Jo 


a  a  r-a  m-a-l+@(n-r+a) 

X  ,C^r-a^m"a+0^n‘r+a^  I  (1"t)  (1"t  '  t  dtHl 

(3.3.4) 


for  r  <  m  +  n,  which  completes  the  proof. 
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An  Integration  by  parte  and  some  manipulations  with  the  binomial 
coefficients  if  required  to  show  that  PQ(a;  m,  n,  r)  as  given  by 
(3.2 .46)  satisfies  the  partial  difference  equation  (3*2.1).  An  example 
of  the  type  of  computations  necessary  to  show  that  the  boundary  condi  - 
tions  are  satisfied  is  given  below.  The  others  require  similar  computa¬ 
tions. 

P0(a;  0,  n,  r) 

/Own.  ,  w  \i  f1  /,  ..a,,  flxr-a  -a.-l+0(n-r+a) 

*  V(r)  I-a+$(n-r+a)]  J  (1-t)  (1-t  )  t  dt 


if  a  >  0 


0  [0(n-r)J 


r 


(i-tTt 


ejr  8(n-r)-l 


dt  if  a=0 


a  >  0 


.(”£T)T  P  (r+l,n-r)  a=0 


a  >  0 


a  =  0 


3.4  Moments:  One  can  always  write  down  an  expression  for  the  kth  moment 
as  a  finite  sum  of  terms  simply  by  using  the  moment  definition.  For 
example,  the  first  moment  is,  by  definition. 


Eg (A;  m,  n,  r)  s 

r  fl  a  a  r-a  m-a-l+8(n-r+a) 

I  aO(r-a)  [n-a+®(n-r-t*)  J  I  (1-t)  (1-t  )  t  dt. 

a"°  ’  (3.4.1) 
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This  is  not  an  attractive  form  and  requires  considerable  computational 
effort  to  obtain  a  numerical  result  for  given  values  of  m,n,r  and  0. 

A  closed  form  of  (3.4.1)  would  be  desirable  but  unfortunately  is  not 
obtainable.  However,  some  rather  interesting  side  results  were  obtained 
as  byproducts  of  attempts  to  simplify  (3.4.1). 

The  first  approach  used  was  to  condition  Eg(A;  m,  n,  r)  on  the 
result  X(l)  of  the  first  draw. 

Eq{Aj  m,  n,  r)  =  E  [eq(Aj  m,  n,  rjx(l)}] 

-[1+VA>  m-1>  n>  r‘1!]  sh  *  ElAf  ">  n-1'  r'11  SSe 

or,  simplifying, 

(m  +  n0)  E0(A;m,n,r)  =  m  Eg{Ajm-l,r-l)  +  n0Eg  {A;m,n-l,r-l]  +  m. 

(3-4.2) 

Clearly,  E_(A;  m,  n,  r)  must  satisfy  the  boundary  conditions 

U 

Eg ( A;  0,  n,  r)  =0  (3*4.3) 

Eg (Aj  m,  0,  r)  =  r  .  (3-4.4) 

Define  D0(m  +  n-r,  m,  n)  =  E0{Aj  m,  n,  r) 


The  (3-4.2)  through  (3-4.4)  become 

(m+ne)Dg(m+n-r,m,n)=mDg(m+n-r,m-l,n)+n0Dg(m+n-r,m,n-l)  +  m 


(3-4-5) 

Dg(n-r,0,n)  =  0 

(3-4-6) 

Dg(m-r,m,0)  =  r  , 

(3-4-7) 

m+n 

Define  „(x)  =  T  Da(m+n-r,m,n)x  . 

m>n  rtl  9 

(3-4-8) 

2  6 


0,»<x>  ■  0 

m 

m 

Wo(x)  *  Tl,  “r  -  * 

I 

r-1 

(3,4.9) 


<L 

d*  V  X-1  / 


mxm+^  -  ( m+1 )  xm+l 


(x-i  y 


(3.4.10 


Multiplying  (3.4.8)  by  (m+n©)  and  using  (3-4.5)  gives 

m+n  m+n  m+n 

(m+n0)£l  (x)  =  £  mD  (m+n-r^m-l^njx  +  T  n0Do(m+n-r,m,n-l)xr+  T  mx 

'  r=l  0  r=l  9  rtl 

Replacing  r  by  r+1  and  noting  the  first  two  terms  are  zero  when  r=0 
gives 

m+n 

<"™e)  awl,nix)  +  “*  nB,n-l<x>  +  “  ' 


(3.4.11 


Nov  define  A(|,tj)  =  T  7  «  •  (3-4.12 

m=0  n=0  m'n  m* 

Then  (3-4.9)  and  (3-4.10)  become 

A  (0,t))  =  o  (3-4.13 

A  (5,0)  ■  L  “%o(x)  £ 


»  _ m+l  ,  ..  i  m  .  N .m 

r-  mx  -I m+l  lx  +1)8 


Differentiating  A  with  respect  to  |  and  r)  and  adding  gives 


I  A.d*1))  +  T)$AT)({,r))  -  £  l  (m+nS)  0  (x)  JA- 

5  m»0  n-0  * 


—  —  ,m_n  oo  oe  /  V  ,iu  u 

T  l  rn  X  a  .  (x)  VV  +  I  I  nexfi  (x'  VL 

A  “-1*13  m!nJ  A  A  “b0’1  mln! 


j  n 


m=0  n*0 


m=0  n=0 


00  00 

I  I 

m=0  n=0 


m  x 


m+n  , 

x _ -1 

x-1 


„m  n 

iA  . 


(3.4.15) 


Taking  the  terms  of  the  right  hand  side  of  (3*4.15)  one  at  a  time 

„m_n 

Y  mxD  ,  (x)  • 

_  Lm.  ▼»>_  i  n'  '  i 

m=0 


oo  oo  m  n  oo  oo  m-1 

r  r  mxn  ,  (x)  ■* ■  §x  y  y  n  (x)  4 — 

M)  n*0  m-:l'n  m!nl  mil  n*)  m'1'n  t53'1 


“-1 

^r,=  xiAd^) 
Tl* 


(3.4.16) 


.m  n-1 


So  S  **  ^  =  9X1  io  £ 

(3.4.17) 

l  y  mx  ^  r  r  (xm+n  .  i^/CX^>-e^ 

mtonto  -1  ^mtlnV  *-<*  / 

(5.4.18) 

Substituting  (5.4.16)  through  (5*4.18)  into  (5*4.15)  and  simplify¬ 
ing  ve  obtain 

e(Aj  -  XA)  +  0tj(A^  -  XA)  =  A  (Cx(5+^-ee^)  .  (5.4.19) 

Let  A(|,ti)  =  ex^+11^  V(|,tj).  Then  (5-4.19),  (5.4.15)  and  (5.4.l4)become 

5  Vg(5,n)  +  env^S/T))  =  ~x  (*  -  e(1'x)(5+^)  (5.4.20) 


V(0,T))  =  0 

v(!,0)  =  —A  bc2 * * 5|-x5-l4€ 
(x-l)2V 


5(1 


-xi>. 


(3.4.2D 

(3.4.22) 
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Let  us  apply  the  method  of  characteristics  to  solve  (3.4.20) : 


at  dn  dV 


e 

The  first  and  second  members  yield  ^  =  c^  . 
The  first  and  third  members  yield 


dV  _  x 

f/v  (l-x)(t+C  t6)')  so  that 

dt  ‘  x-1 

^x-e  1  j 

0 

II 

-P 

-P 

£dt  =  dx.  Then 

vU,») 

(i- .( !-*)(«»  ♦»**)„. 

'  o 

We  see  that  the  boundary  conditions  (3. 4.21)  and  (3-4.22) 
since  V(0,p)  =  0  and 


v(!’°>  ■  -TTT5JT  ['(1'X)5-1]) 

x  f  2.  ,  .  (l-x)t 

- - 5  x  |  -  x£  +  ev  -lJ 

(x-1)2  V  / 

Therefore 

A(t,n)  =  -  f  e(l-x)(T|t  +  ^it) 

\  •'o 

=  ~l  Z  X  (xm+n+1-f  [(l-x)t+x]m[(l-x)tl 
x  x  m=0  n=0 V  J0 


(3-4.23) 


(3-4.24) 
are  satisfied 


Which  yields  flm^n(x)  as  the  coefficient  of  ; 


m  n 

A 

min: 


flm,n<x>  *  Sr  (*+n-  f1  [(l-x)t+xr1[(l-x)t^]ndt) 

\  /(3.4. 


(3.4.26) 


To  recover  D@(m  +  n  -  r,  m,  n)  we  note 

r*  D0(m+n-r,m,n)  =  (a  (x)\  .  Prom  (3-4.26) 

dx  \  '  /  x=0 


and 


r 

E 

J=o 


-  E  (J) 


fer 


I  (l~t)x+t  j  [(l-te)x+t0]  dt  \  ( ~ 

dxr"^  '1_x> 


Now 


dr-J 


dxr'J  \1-x/x=0 


fer 


0  J  -  r 

(r-j)!  J  <  r 


lx=0 
(5.4.27) 

(3-4.28) 


m_1  a  a  n 

[(l-t)x+t]  [(l-te)x+tS]  dt 


.] 

Jx=0 


-  I  <J> 

k=0  * 


I1  fe 

Jo  |_dx 


f(l-t)x+t]  — r-^  [(l-te)x+t6]  Idt 

dxJ"  I 

Jx=0 


-  jo  (J)(«-i)k(„)J_][  J1  t”-1-kti.t)kte(n-',+k>U-te)l£-J 


dt 


=  t  J!  (mk1)(j?k)  f\l-t)  (l-te)k->Jtm-1"k+«(n-J+k)( 

k=0  -'o 


dt 

(5.4.29) 
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Substituting  (3.4.28)  and  (3.4.29)  into  (3-4.27)  gives: 


,(Aj»,n,r)  .  m  rf  £  (“t1)(A)  f  ( 1-t)  (1-t9) 

J-0  k«0  Jo 


1  k  fl  J-k  m-k-l+0(n-J+k) 
t  dt 


(3.4.29) 


which  is  certainly  no  improvement  over  (3.4.1), 


Another  approach  is  to  recognize  that 

1  dr 


E  (A>  m,  n,  r)  -  — 

e  r.  ^r 


y=lJx=0 


where  4>m  n  is  defined  by  (3.2.19).  Using  (3.2.43)  we  find 


&3  (x,y) 

m.n'  ’  ' 


~Qy" 


=  Hi 


ly=l 


m-l)  f 1 

r'.  J 

Jo 


m-2 


0 


((l-x)t  +  x]  [(l-x)t  +  x]  x(l-t)dt 


~  f  [(l-x)t  +  x]  [(l-x)t0+  x]  x(l-t)t0~1dt 

r-  Jo 

(3.4.30) 

Denote  the  first  term  on  the  right  hand  side  of  (3.4.30)  by  h(x)  and 
apply  Leibniz's  rule: 


arh(x) 


dx 


3  Hte^l  r1  <£_  x  [(1.x)t  +  x]m  2  [(l-x)t0  +  x)  (l-t)dt 

.  ‘  Jo  dx 

x=0 

■  ■ ■%“  r  i®  i s-'-rf-wo-'L 


f  r  ^IT  [(I"*)1  +  x]Jm2^(l-x)t0+  xJV-t)dt 
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■ «  /.  I  «‘>  !?['->■  •  *r  S^"'*  • 

.  iVs1)  J'l(^ikt^-,‘(i-t)l"1(n)r.1./(n-"1’k)i-tVl-‘« 


r-1  _  _  /I  k+l  9 

m(m-l)  X  (“k  )(r-8-l)  (1-t)  (1-t  ) 

k=0  Jo 


1  k+1  9  r-k-1  m-2-k+0(n-r+k+l) 


t  dt . 

(3.4.31) 

Similarly,  the  second  term  on  the  right  hand  side  of  (3.4.30), 
when  differentiated  r  times  and  x  placed  equal  to  zero  gives 


r-1  /l  k+1  6  r-k-1  m-2-k4£(n-r+k+l) 

l  (VXrStil)  I  (1-t)  (1-t)  t  dt. 

k=0  •Jq 


mn0 


Combining  (3.4.31)  and  (3-4.32)  gives 

E  (A;m,n,r)  =m  X  i (m-1) (™kZ)  (r-ic-i) 

9  k=0  \ 


(3-4.32) 


.  ‘j  (1  k+1  9  r-k-1  m-2  -k+0  ( n-r+k+1 ) 

“■1\  /  \!  I  /1  x\  -/l  X  \  J L  ^ 

(3.4.33) 


j  (l  k+1  8 

+  n0(\1)(r^kil)|J  (1-t)  (1-t)  t  dt. 


k+1  k  k 

Writing  (1-t)  =  (1-t)  -  t(l-t)  and  breaking  up  the  integral  in 

(3.4.33)  into  two  parts  we  obtain 


E0(A;  m,  n,  r) 


r-1  -1  n  „  n-1  f1  4  q  r-l-k  1 
-X  k{m(mk1)(r-l-k)+  n0(S)( r-l-k)}  (1-t)  (1-t  )  t 
k=0  'Jo 


n-1  fl  k  q  r-l-k  m-k-l+0(n-r+l-k) 

dt 


r-1  n  fl  k  9  r-l-k  : 

+  l  (V) (r-l-k)  (1-t)  (1-t  )  t 

k=0  Jo 


k  9  r-l-k  m-k-l+0(n-r+l-k) 


dt 
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E@(A;  m,  n,  r-1) 


+  m 


*1  ("^Kr-i-k)  f1  (i-t)VtV'1'^"'k‘l48("‘r+1'ltV 

kaO  J0 


(3-4.34) 


r 

But,  since  A  =  ][  X  (j),  we  have 
J=1 


Eg{A;  m,  n,  r]  =  E0(A;  m,  n,  r-1}  +  P(X(r)  =  1)  and  hence,  in  (5 *5^) > 
we  can  identify  the  unconditional  probability  of  drawing  an  element  from 
M  on  the  r  draw  as 


P{X(r)=l} 


m  ^  (Vl(r-i-k) 

k=0 


(l.t)k(l-te)r-1“ktm-k-1+e(n'r+1-k)dt 

(3.4.35) 


which  is  a  valuable  result  itself. 

Equation  (3-4.35)  gives  us  insight  into  interpreting  (3-4.29) 
Replacing  j  by  j-1  in  (3.4.29)  gives 


EQ(Aim,n,r)  =  m 


j=l  k=0  Jo 


r 


=  I  P(X(J)  =  1) 

J«1 


so  that  (3.4.29)  simply  states 


r 

E  (A;  m,  n,  r)  =  £  E  [X( J) ] . 

8  jti 


SECTION  4 


ASYMPTOTIC  RESULTS 


4.1  When  m  and  n  are  large  compared  to  r. 

Using  (3*1.4)we  shall  obtain  crude  bounds  on  Pg(A;  m,  n,  r)  which 
converge  to  the  same  quantity  as  m  and  n  go  to  infinity  while 
m/n  =  ti  is  held  fixed.  We  first  prove  a  lemma  concerning  each  term  on 
the  right  hand  side  of  ( 3*1.4) v 

Lemma  4.1.1  For  all  k  and  J  such  that 

1  <  k  <  r  and  1  <  J  <  (a)>  6  1  implies 


r  <  r< 

m(k,  x^r^)+0n(k,  x^)  “  m(k,Xj)  +  0n(k,Xj)  "  m(k,x1)  +  0n(k,x1)  § 


Proof:  Let  0  =  1+  c,  €  >  -  1.  For  all  k  and  j  such  that 
l<k<r,  1  <  J  <  (a) 

m(k,Xj)  +  0n(k,Xj)  =  m(k,Xj)  +  n(k,Xj)  +  en(k,Xj) 

=  m+n-k+1  +  €n(k,Xj)  by  (3.1.2)  . 


Sy  the  definition  of 


n(k,Xj) 


and  the  ordering  imposed  on  x(a)  we  have 


n(k,  x(g))  5  n(4,Xj)  <  nU,*^  . 


Multiplying  this  inequality  through  by  e  (the  inequalities  are 
reversed  if  €  <  0)  and  adding  m+n-k+1  to  each  term  proves  the  lemma. 


Let  n  =  m  +  n0  and  (d)  =  d(d-l) . . .(d-m+l) . 

m 

Using  (3.1.2)  and  the  definition  of  x1  we  see  that 
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TT  m(k,x1)+0n(k,x1)J  =»£n(^-l) . .  .(n-a)(n-a-0) . .  .(n-a-[r-a-l]0)j 


» 


<*>a  <¥>r-a  •M 


-1 


(4.1.1) 


<s> 


and,  using  the  definition  of  x^ 

77  m(k,x^^)  ©n(lt,3C(^))^  =  . .  .(n-[r-a]0)(n-[r-a]0-l) . . . 

(n-[r-a]0-a+l) j 


•  [<$> 


_  (ti-[r-a]0)  e 

r-  a  a 


-a]'1 


(4.1.2) 


Using  lemma  4.1.1  and  substituting  (4.1.1)  and  (4.1.2)  into  (3.1.4)  for 
0  1  gives 


_  (m)  (n) 

(r)  a  r~a 
a  (m+n0)  (n+  ^3-) 


0  r-a 


Pg{a;m,  n,r) 


(!) 


(m)_  (n) 


a  r-a 


a  (m+[n-r+a]0)&(n+  ^)r_a 


(4.1.3) 


For  0=1,  the  bounds  coincide  and  equal  the  correct  central 
hypergeometric  probability. 

Lemma  4.1.2 

(am  +  b) 

lim  7 - -  -rr  =  (-)  where  a,b  and  d  are  real,  a  is  any 

(cm  +  d)^  c  ' 

m-»oo  'a 

positive  integer  and  c  ^  0. 

Proof:  The  proof  is  immediate  by  writing: 
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lim  (“+b)a 


11m 


b-a+i 


m 

m-».  (cm+  d)a  -  m—  (c  +  £)(c  +  ^)(c  +  _(c  +  dHDj+lx  =  * 

mm  m  m 


) 


Let  1  =  “  in  (4.1.3)  and  apply  lemma  4.1.2  to  both  the  upper 


bound  U  and  the  lower  bound  Lj 


(m)  (rim) 

U  =  (T)  - 2 - ^ - - 

([1+Ti0]m  +  [a-r]0)a(h+  ^m^ 


->  ^  ^1-kj0^  ^1+ne^  (4.1.4) 


(m)  (nm) 

L  =  (T)  _ &  r~a _ 

a  ( [1  +  *l0]m)(  [■£  +  Tj]m  -  -|) 


_>  (a}  ClTn^)  (l^e)  (4*1*5) 


r-a 


so  that  the  bounds,  and  hence  Pg( a;  m,  n,  r),  both  approach  a  binomial 
probability  with  parameters  r  and  p  =  y^'g'  .  Therefore,  if  r  is 
small  with  respect  to  m  and  n,  one  can  use  the  properties  of  the 
binomial  distribution  to  investigate  pertinent  questions  such  as  esti¬ 
mation  and  moment  problems.  For  example,  the  uBual  minimum  variance 
unbiased  estimate  of  p  in  a  binomial  distribution  can  be  used  to 
estimate  9  by 

?  =  P  =  l^I  3=  f/y—j  (4.1.6) 

and  the  first  two  moments  of  the  distribution  P0( A;  m,  n,  r)  are 


rm 

11  m  +  n0 

2  r  m  n  0 

o  =  - p 

(m  +  n0) 


(4.1.7) 


For  m  =  n  =  20,  table  4.1  shows  the  comparative  values  of  the 
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exact  probabilities  Pg  (a,  20,  20,  r)  with  the  binomial  approximation 
B(p,  r)  with  p  =  . 


Table  4 . 1 


B 

a 

P0(a;  20,20,r) 

B<i^> 

r) 

CVJ 

11 

Q> 

0=4 

0=2 

0  =  4 

3 

0 

1 

2 

3 

.2808 

.4599 

.2260 

.0333 

■ 

1 

.2963 

.4444 

.2222 

.0370 

.5120 

.3840 

.0960 

.0080 

0 

.1088 

.2918 

.1317 

•  3277 

1 

.3278 

.4268 

•  3292 

.4096 

5 

2 

•  3534 

.2234 

.3292 

.2048 

3 

.1704 

.0523 

.1646 

.0512 

4 

.0368 

.0055 

.0412 

.0064 

5 

.0028 

.0002 

.0041 

.0003 

4.2.  When  m,  n,  and  r  are  large. 

Let  us  now  investigate  the  behavior  of  PQ[A;  m,  n,  r)  as  given  in 

(3.2.46)  as  m,  n  and  r  approach  infinity  while  the  ratios  “  = 

and  —  =  p  are  held  fixed, 
m 

Lemma  4.2.1  (Laplace's  method) 

Let  h(t)  be  a  twice  continuously  differentiable  real  function  of  the 
real  variable  t  on  a  <  t  <  P  and  let  g(t)  be  continuous  on 
a  <  t  <  P.  Let  h(t)  have  a  unique  maximum  at  t  =  tQ,  0!  <  t  <  P, 
and  g(tQ)  ^  0.  Then,  for  large  m, 

/e(t)e"(t)dt~/^  g(t0)e  °}  (4.2.1) 

a 
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Remark:  0  If  h{t)  has  a  finite  number  of  maxima,  the  integral  can 
be  broken  into  a  finite  sum  of  integrals,  each  satisfying  the  unique 
maximum  requirement. 


d)  If  tQ  =  P  and  h(t)  is  well  behaved  at  t  =  P,  then  the 
right  hand  side  of  (4.2.1)  must  be  halved. 

Lemma  4.2.2 


Let  ^  =  p  >  0,  —  =  T|  >  o  and  t  = 

m  m  o 


1+011 


1 

e 


Then,  for  large  m. 


A 

J  (i-t)a(i-tV"atm'a_1  +  0(n_rfa)dt 


(l-tc)  (1-t-) 


r-a+1  m-a+0(n-r+a-  -|) 


Proof:  Assume  the  parameters  have  been  labeled  so  that  r  <  n.  Thus  can 
always  be  done  without  loss  of  generality  and  insures  0  <  tQ  <  1. 
Rewriting  the  integral  in  the  form 


a  a(0-l)-l 

where  g(t)  =  (  ~\)  t 
1-t 


Q 

and  h(t)  =  pin  (l-t  )  +  [l+0(vp)  ]in  t,  we  note  that  h(t)  is  twice 
continuously  differentiable  on  (0,l)  and  g(t)  is  continuous  on  (0,l). 
Furthermore,  h'(t)  =  0  has  the  unique  solution 


i+e  (n-p) 


1 

10 


o  =  [  1+0T1 


p02t0 

and  h"(tQ)  =  -  — ^ - 


tc(l+t0)2 
o'  o' 


<0  so  that  t  is  the 
o 


unique  maximum  of  h(t)  on  (0,l).  Lemma  4.2.1  therefore  applies  and, 
after  some  simplification,  the  stated  result  is  obtained. 


Lemma  4.2.3  For  —  =  p  >  0,  t  1  >  0 
1  1  — m  m 
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Proof: 


(i) 

(ii) 


a  - 

(m-a) 


1 

2 


n 

n 


+ 


1 

2 


n-  r+a+ 

VS  (m-r+a) 


^(r-a) 


n-a+ 


1 

2 


( i)  By  Stirling's  formula 


1_ 

a! 


(5) 

_______  m-a 

V2*(m-a)(—^) 


/m_  e-a 
m-a 


-m 


(l  -  ^)  (m-a) 

m 


a 


_1 

a! 


1 

m2  (m-a) 


1 

2 


(ii) 


/  n  ■>  T)ml 

^r-a'  “  (pm-a) I(  Fl-p]m+a) ! 


; -  (H2) 

/2nT)m  e' 


^m 


Fn-pTm+a 


_  pm-a 

V2n  (pm- a)  (£5"-)  /2fl(  [n-p]m+a) 


( [yp]m+a^ 


_  run  /  \  pm  a 

1  /  ~n  /  nm  x  /(T]-p)m+a\  /  pm-a  x 

^  ( r-a)  ( n-r+a)  '  (T[-p)m+a'  v  pm-a  '[ii-p)m+ay 


,  _ _  n  .  r  a 

1  n  /  n  \  /n-r+as  ,  r-a  \ 

^  ( r-a) (n-r+a)  'n-r+a'  r-a  '  'n-r+a' 


n-r+a+  4  r-a+ 
(n-r+a)  (r-a)  ^ 
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Now,  applying  lemmas  4.2.2  and  4 .2 . 3>  we  obtain  the  asymptotic  expansion 


Pfl  (a;  m,  n,  r)  ~ 


|m-a+0[n-r+a]\  1 

1  nm  t® 

(  0( n-r+a]  I 

r(r-a)(m-a)(n-r+a) 

(4.2.2) 


8iV  ...w.  .8  -r+1 


X 


/( m-a)(r-a)(l-to)t^  ^i-r+a)  (l-t®) 
^n-r+a)(l-t®)to 


fi-r+a)  (] 
,  (r-a)t^ 


0  \n 

nto|  tm 

- — !  o 


n-r+a 


4.J.  When  m,  n,  r  and  a  are  large : 

Using  lemma  4.2.1  we  now  prove 
Lemma  4.3.1. 

j*  n  fl 

Let  —  =  p  >  0,  -  =  T|  >  o,  —  =  a  >  o.  Then  there  exists  a  unique 

solution  t^  e  (0,l)  of  the  equation 
0+1  0 

(1+  q)t  -  (l+@T)-a)t  -  (l+0[ri-p-K2])t+l-a+e[q-p+a)  =  0 


and  for  large  m. 


I 


1 

(i-t)a(i-te)1”atm-a-1+e(n-r+a) 


dt 


2ir 


\ 


at1(l-te)2+(r-a)02t®(l-t1)2 


( l-tx)  a+1(  !-t®  )r-^ltlm-a+0(  ^r+a) 


(4.3.1) 


Proof:  Rewriting  the  integral  in  the  form 

f  where  g(t)  =  f1  and  h(t)  -  /h(l-t)a(l-te)p-at1-0M[r|-pW:1> 

J  0 

we  note  that  g(t)  is  continuous  in  (0,1)  and  h(t)  is  twice  continuosly 
differentiable.  Furthermore,  h(o)  =  h(l)  =  -  »  so  that  there  exists 
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at  least  one  relative  nmviimim  off  h(t)  in  (0, l) .  Let  t^  be  any  solution 
in  (0,1)  of 

h.(t) _ SL  .  (p-a)et9;1  +  i.-0f^[T)-P^  =  0_ 

1"t  l-t6  1 

Simplifying,  let  t^  be  any  solution  of 

(i+6T))t0+1  -  (i+0n-a)t9  -  (i+e[T)-p+a])t+i-a+e[n-p+a]  =  o,  (4.5.2) 


in  (0,l).  To  show  that  there  is  only  one  such  solution,  it  suffices  to 
show  that  h”(t^)  <  0  since  a  continuously  differentiable  function 
cannot  have  two  relative  maxima  without  a  relative  minimum  separating 
them. 


h-'(t)  =  .  a  .  (p-a)et9~2(e-i+t9)  _  i-a+e(yp-Kx) 

(l-t)2  (l-t6)2  t2 


h,,(t1)  = 


a 


(p-a)02t9-2 

(l-t9)2 


(p-a)et 


e-i 


l-t 


e 


1 _ l-a+ed-p-ta) 


t. 


a 


(l-tl)< 


(o-a)62t®'2 

(l-t9/ 


a 


TT^ 


a  (p-a)02t9 

d-ti)2ti  d-t9)2t2 


since  p  -  a  >  0 


and  0  <  t^  <  1.  Hence,  t^  is  the  unique  maximum  of  h(t).  Applying 
lemma  4.2.1  then  proves  the  lemma. 

Using  Stirling's  formula  we  obtain 


m 

2na(m-a) 


(-2L) 

m-a 


4l 


m 


a 


(4.5.5) 


Tables  4. 3-1  and  4-3-2  furnish  examples  of  the  accuracy  of  the 
asymptotic  expansion  (4-3-5)  for  moderate  and  large  values  of  the 
parameters,  respectively. 
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Table  4.3.1 
m  *  n  *  20,  r  =  14 


e 

-  2 

e 

-  3 

Exact 

Approx. 

A 

Exact 

Approx. 

A 

mm 

— 

— — 

.00245 

| 

.00041 

■02579 

.02812 

.03040 

.03186 

.00146 

.10502 

•10999 

.10478 

.10832 

.00354 

.22253 

•22993 

.00740 

.21141 

•21713 

.00572 

.27664 

.28400 

.00736 

.26587 

.27211 

.00624 

. 21483 

.21981 

.00498 

•21585 

.22054 

.00469 

•10793 

.11025 

.00232 

.11521 

•11767 

.00246 

•03572 

.03647 

.00075 

.04066 

.04156 

.00090 

.00783 

.00800 

.00017 

.00943 

.00966 

.00023 

.00113 

.00116 

.00003 

.00141 

.00145 

.00004 

.00011 

.00011 

0 

.00012 

.00014 

.00002 

0 

0 

0 

0 

.00001 

.00001 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Table  k-3-2 

m  *  95,  n  =  120,  r  =  75 


SECTION  5 


EXAMPLES 


5*1  Example  1:  A  Genetics  Problem 

Let  us  analyze  the  data  in  table  2.2.  Using  the  normal  approximation  to 
the  central  hypergeometric  distribution  with  parameters  m  =  95*  n  =  120, 
and  r  =  75*  Fisher's  exact  test  requires  rejection  of  the  hypothesis 
of  unbiased  selection  at  the  5^  (nominal)  level  if,  and  only  if. 


>  K 


.025 


rmn(m+n-r) 
k(m+n)2(m+n-l)j 


1 

2 


where  K  =  1.96  is  the  97-5  percentile  of  the  unit  normal  distri-.' 
bution.  For  the  values  appearing  in  table  2.2  we  find  la  -  =7.9 

1  w+n  1 

and 

K  I  rmn(mtn-r) 

,025  \(m+n)2(m+n-l) 

so  that  the  null  hypothesis  is  rejected.  If  we  define 
4° 

f(©)  =  ^  Pg(a;  95,  120,  75) ,  we  can  display  the  sensitivity  of  the  above 
a=27 

test  in  the  form  of  an  operating  characteristic  curve.  For  example, 
referring  to  Figure  5.1.1*  we  see  the  probability  of  accepting  the 
hypothesis  of  randomness  when  indeed  "nature"  is  "twice  as  hard"  on  the 
heterozygotic  rabbits  (i.e.,  ©  =  2)  is  0.17. 


PROBABILITY  OF  ACCEPTING 


1.0 


Suppose  a  lot  of  transistors  is  to  be  divided  among  two  purchasers.  The 
lot  is  dichotomized  by  some  quality  criterion  and  let  us  assume  thfere 

k6 


are  m  high  quality  items  and  n  other  items  of  tolerable  quality. 
Purchaser  G,  who  has  contracted  for  r  items,  subjects  the  supplier 
of  favoring  purchaser  H  with  respect  to  the  distribution  of  the  m 
high  quality  items.  He  therefore  decides  to  test  the  hypothesis  that  the 
lot  has  been  divided  in  a  random  manner  (definition  l.l)  against  the 
one-sided  alternative  that  the  supplier  is  biased  (in  the  sense  of 
definition  2.1)  in  favor  of  purchaser  H.  Fisher's  exact  test  is 
appropriate  and  the  number  of  high  quality  items  received  by  purchaser 
G  has  a  noncentral  hypergeometric  distribution  with  parameters  m,  n, 
and  r.  In  terms  if  the  noncentrality  parameter  0,  purchaser  G 
wishes  to  test 

H  :  0  <  1 

o  — 

:  0  >  1 

Let  A  be  the  number  of  high  quality  items  received  by  purchaser  G  in 
a  shipment  of  r  items.  Hq  is  then  rejected  if,  and  only  if,  A  <  Aq 
where  Aq  depends  on  m,  n,  r  and  the  significance  level  a.  Figures 
5.2.1  -  5-2.5  show  operating  characteristics  for  m  =  n  =  10  (a  »  .10), 
m  =  n  =  12  (a  «  .11)  and  m  =  n  =  20  (a  «*  .05) •  The  curves  are  indexed 
by  various  values  or  r  and  the  required  Aq. 

5.5  The  Noncentral  Hypergeometric  Distribution  and  a  Class  of  Urn  Problems. 
In  this  subsection  we  shall  consider  a  class  of  urn  sampling  problems 
which  lead  to  the  noncentral  hypergeometric  distribution.  Suppose  an  urn 
contains  m  type  1  chips  and  n  type  2  chips.  A  chip  is  drawn  at 
random  and  kept  with  probability  p^  if  it  is  type  i,  i  =  1,2.  With 
probability  1  -  p^^  the  drawn  chip  is  returned  to  the  urn  and  a  second 
chip  is  then  taken  at  random.  This  process  is  repeated  until  a  chip  is 
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Figure  5*2.1. 
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Figure  5 • 2 . 3 • 


retained.  Let  X  and  Y  be  random  variables  defined  by 


X  =  i  if  the  first  drawn  is  a  type  i  chip 

Y  =  i  if  the  chip  retained  is  a  type  i  chip. 

i  =  1 


Clearly 


P{X  =  i)  = 


m 

nH-n 


n 

m+n 


1=2 

otherwise 


and  P{Y=l}  is  the  probability  that  a  type  1  chip  is  drawn  and  retained 
before  a  type  2  chip  is  drawn  and  retained.  Hence, 

P{Y=l}  =  E[P{Y=1|X}] 

=  P[Y=1 |X=l)P{X=l)  +  P{Y=1 |X=2 }P{X=2 ) 

-  (Pl  *  <i-Pl>m.i)i  ^  +  [d-p  )p(i-D]  ^ 


Solving,  we  find 


P{Y*1)  = 


m 


and  observe  that  the 


m+(~)n 


*2 


sampling  mechanism  has  bias  6  =  —  in  the  sense  of  definition  2.1. 

Therefore,  if  one  continues  sampling  in  manner  described  above  until 

r  chips  have  been  retained,  the  number  A  of  type  1  chips  in  the  sample 

has  a  noncentral  hypergeometric  distribution  with  parameters  m,  n,  and 

r  and  noncentrality  depends  only  on  the  ratio  of  p2  to  p^^  so  that 

the  distribution  of  A  is  unchanged  if,  for  example,  p^  >  p2  and 

type  1  chips  are  kept  with  probability  1  while  type  2  chips  are  retained 
*  p2 

with  probability  p  =  — . 
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U.  S.  Naval  Boiler  and  Turbine  Lab. 
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Commanding  Offlcar  and  Director 
U.  S.  Navy  Electronics  Lab. 
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U.  S.  Naval  Propellant  Plant  <QS) 

Indian  Head ,  Maryland  1 


Commander 

Southern  Air  Materiel  Area,  Pacific 
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Mr.  G.  Burrows 
Knolls  Atomic  Power  Lab. 
6  Schenectady,  New  York 


Technical  Operation,  Inc. 

andC.  0.  R.G. 

Hq.  Continental  Army  Command 
Fort  Monroe,  Virginia 


Dr.  A.  Charms 
The  Technological  Institute 
Northwestern  University 
1  Evanston,  Illinois 


U.  S.  Naval  Avionics  Facility 
Indianapolis  18,  Indiana 
Attn:  Library 


U.  S.  Naval  Engineering  Experiment 
Station 


Anapdis,  Maryland 

Attn:  Mr.  F.  R.  DelPriare 


Mr.  W.  H.  Clatworthy 
Bettis  Plan,  Wesllnghouso 
Electric  Corporation 
Box  1468 
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